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Abstract  -  A  novel  imaging  mode  of  cone-beam  volume  CT  is 
proposed  in  this  paper.  It  adopts  a  raster  scanning  x-ray  source 
on  a  trough-like  surface,  and  a  group  of  point  detectors 
distributing  on  a  large  circle  plus  an  orthogonal  arc.  Through  a 
single  rotation  of  the  trough-like  surface,  a  full  set  of  projection 
data  can  be  acquired  and  an  accurate  3D  image  can  be 
reconstructed.  This  paper  describes  the  shape  and  dimension  of 
the  trough-like  x-ray  source  and  the  configuration  of  detectors 
in  detail,  proves  that  this  structure  is  satisfied  with  the 
completeness  condition,  and  gives  a  reconstruction  algorithm  of 
3-D  image  adapting  to  this  structure.  Computer  simulations 
show  that  this  structure  could  achieve  optimal  result  in  meeting 
the  following  requirements:  minimum  time  of  rotation,  less 
Compton  scatter,  and  fast  image  reconstruction  speed  while 
keeping  the  completeness  condition  to  be  realized. 

Keywords  -  Volume  CT,  completeness  condition,  trough-like 
surface,  point  detector,  image  reconstruction 

I.  Introduction 

As  an  extension  technology  of  traditional  computed 
tomography,  direct  volume  CT  has  been  a  hot  topic  in  recent 
years.  Comparing  with  traditional  CT,  volume  CT  improves 
the  dose  usage  of  X-ray,  the  speed  of  projections  collecting, 
and  most  important  it  is  able  to  provide  isotropic 
reconstruction  images  in  spatial  domain. 

To  reconstruct  a  3-D  image  exactly  from  2-D  projection 
data,  Tuy  [1]  showed  that  each  projection  plane  passing 
through  the  object  should  intersect  the  orbit  of  the  source 
vertex.  This  requirement  is  referred  to  as  Tuy’s  data 
completeness  condition.  Only  if  the  completeness  condition 
is  met,  the  precise  reconstruction  could  be  ensured. 

At  present,  there  are  several  methods  for  direct  volume 
CT,  such  as  direct  volume  CT  based  on  ideal  vertexes  of  the 
cone-beam  [2],  direct  volume  CT  based  on  limited  projection 


[3],  helical  or  spiral  volume  CT  [4],  These  scanning  modes  of 
cone-beam  CT  have  the  following  shortcomings:  (l)the  ideal 
non-planar  geometry  of  the  vertex  orbits  are  almost 
impossible  to  be  implemented,  (2)the  effect  of  Compton 
scatter  to  planar  detectors  is  so  heavy  that  reconstruction 
images  have  severe  artifacts,  (3)the  reconstruction  algorithm 
is  very  complex  and  time-consuming,  and  so  on. 

In  this  paper  we  break  all  above  conventional  scanning 
structures  of  volume  CT,  and  propose  a  new  imaging  mode  of 
volume  CT  based  on  a  trough-like  x-ray  source  and  point 
detectors.  This  structure  could  gain  full  projections  through  a 
single  time  rotation  of  the  trough-like  x-ray  source. 

II.  NEW  imaging  mode  of  Volume  CT 

A.  Reverse  G  eometry  of  X-ray  Radiology 

R.  D.  Albert,  T.  M  .  Albert  proposed  a  reverse  geometry  of  x- 
ray  (RGX )  technology  [5].  The  principle  of  the  RGX  system  is  that 
the  object  is  placed  next  to  a  large  raster  scanning  x-ray  source  with 
a  point  detector  located  5cm~3m  away  from  the  object. 

The  remoteness  and  small  size  of  the  detector  eliminates 
nearly  all  scattered  x-rays  from  detection,  therefore  it  would 
improve  resolution  and  avoid  artifacts  in  the  reconstructed  images. 

B.  A  Direct  Volume  CT  with  X-ray  from  a  Trough-like  Surface 
and  P  oint-detectors  on  C  ircl e-plus-arc  C  urve 

We  adopt  RGX  mode  in  direct  volume  CT  and  mend  the 
structure  of  RGX  in  following  ways.  Its  scanning  structure  is  shown 
in  Figure  1.  Firstly,  to  gain  complete  projections,  we  use  a  trough¬ 
like  surface  instead  of  a  planar  raster  scanning  x-ray  source. 
Secondly,  in  stead  of  a  single  point  detector,  we  use  a  group  of  point 
detectors  distributed  on  a  large  circle  plus  a  perpendicular  arc  in  the 
gantry  to  collect  projections,. 
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Figure  1.  The  structure  of  a  novel  volume  CT  with  x-ray  from  a  trough-like  surface  and  point-detectors  on  circle-plus-arc  curve. 
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The  point  detectors  which  evenly  distribute  on  a  large  circle 
and  an  arc  are  indicated  by  open  circles  in  Figure  1.  The  trough-like 
x-ray  source  at  the  initial  position  is  to  be  rotated  to  the  second 
position  as  indicated  by  the  gray  shape  in  the  figure  to  ensure  the 
completeness  of  projections.  To  simplify  the  mathematical 
derivation,  it  is  assumed  that  the  object,  the  x-ray  source  orbit  and 
the  detectors  are  concentric  circles  with  the  center  at  point  0,  and 
radii  r1(  r2,  r3  respectively.  They  should  meet  certain  conditions 
according  to  Tuy's  completeness  condition. 

According  to  the  completeness  condition,  x-ray  from  the 
trough-like  surface  should  not  be  sheltered  by  itself.  Through 
calculating,  we  have  a  set  of  approximate  solutions:  r2=3r1, 

r3=8r1. 

Now  let’s  show  the  length  of  the  perpendicular  arc,  the 
other  part  of  detectors  orbits. 

To  satisfy  the  completeness  condition,  and  in  terms  of 
the  above  set  of  approximate  solution,  we  set  the  radian  of 
perpendicular  arc  to  4arcsin|. 

Furthermore,  we  need  to  know  the  width  of  the  trough-like 
surface.  Through  the  extreme  point  of  the  perpendicular  arc  we  draw 
a  tangent  line  of  the  object.  The  point  where  the  tangent  line 
intersects  the  trough-like  surface  is  the  extreme  point  of  the  trough. 
After  calculating,  we  set  its  width  to  3^ . 

C.  The  Proof  of  the  Completeness  Condition 

Through  a  single  movement  of  the  trough-like  x-ray  source, 
the  volume  imaging  system  described  in  this  paper  could  "wrap"  the 
object  completely.  As  is  known  that  any  plane  through  the  object 
intersects  the  orbits  of  the  cone-beam  x-ray  at  least  once.  So  the 
system  satisfies  Tuy's  completeness  condition. 

M  oreover,  from  the  view  of  the  detectors,  they  evenly 
distribute  on  a  large  circle  and  a  perpendicular  arc,  which  is 
equivalent  to  that  the  vertex  of  a  cone-beam  x-ray  moves  on  the 
orthogonal  circle-and-arc  orbit.  This  shows  that  the  scanning 
structure  satisfies  completeness  condition.  We  have  to  point  out  that 
the  detectors  at  the  lower  part  of  the  large  circle  can't  gain  all 
projections  for  sheltered  by  the  trough-like  surface  itself.  But  the 
radian  of  perpendicular  arc  is  4arcsin|,  then  the  detectors  on  the 

perpendicular  arc  can  gain  the  projections  which  are  sheltered.  From 
all  of  above,  the  direct  volume  CT  system  satisfies  Tuy's 
completeness  condition. 

III.  Image  Reconstruction 
A.  Cone-beam  P rojection  and  3-D  Radon  Transform 

3D  Radon  transform  and  Grangeat's  formula  are  the  typical 
and  efficient  reconstruction  algorithm  for  volume  CT  based  on 
orthogonal  circle-and-arc  orbit. 

The  cone-beam  projection  g(S,A )  of  the  object  function 
fix)  can  be  expressed  as 

g(S,A)=fj(OS+ta)dt  (1) 

where  g(S,A )  is  the  projection  from  the  cone-beam  focal 
point  S  to  a  point  A  in  the  detector  plane. 

The  Radon  transform  of  the  object  is  defined  as  the  plane 
integral  of  f(x) .  In  Figure  2,  any  Radon  plane  C,  can  be  defined 

by  a  unit  vector  0  and  a  scalar  p  where 


Figure  2.  The  Radon  transform  is  the  plane  integral  of  the  object  f(x) . 

<9  =  (sin<2cos<p,sin<9sin<p,cos6>)  (2) 

is  the  normal  vector  of  plane  £( 9,p )  and  p  is  the  distance  from 
this  plane  to  the  origin  of  the  coordinate  system.  The  Radon 
transform  of  the  object  f  (x)  is  given  by 

R(0,p)  =  Hlf(x)S(x-0-p)dx  (3) 

where  the  S function  constrains  the  3D  integral  within  plane 
W,p). 

So  the  object  function  f(x)  can  be  exactly  reconstructed  by 
using  the  inverse  3D  Radon  transform  [6] 


where  dco  =  sin  .  f(x)  can  be  reconstructed,  if  we  could 

gain  R'(0, p)  =  —R[9, p)  from  the  projection  data  g(S,A). 

dp 

Grangeat's  formula  can  complete  this  conversion  process. 

8.  Grangeat's  Formula 

Grangeat's  formula  relates  the  cone-beam  projection  of 
f  (x)  with  the  first  derivative  of  its  3D  Radon  transform. 

In  the  local  coordinate  system  (p,v,a>)  of  Figure  3,  the 
virtual  detector  plane  £  is  defined  by  the  uaxisand  co  axis.  The 

p  axis  is  coincident  with  05  .  The  Radon  plane  £(0,p),  where 
the  plane  integral  takes  place,  goes  through  the  focal  point  S  and 
intersects  the  detector  plane  £  at  lineD1D2.  Another  local 
coordinate  system  [p,  p,q)  is  defined  with  the  rotation  of  the  v 
axis  and  co  axis  about  the  /J  axis  by  an  anglea  (ae[0,;r)), 


Figure  3.  The  coordinate  system  of  cone-beam  projection  to 
establish  Grangeat's  formula. 


such  that  the  p  axis  is  perpendicular  to  line  DjD2  .  The  distance 
from  the  origin  of  the  coordinates  to  line  D  XD  2  is  p.LetpointA 
be  a  point  on  line  DjD2  with  (0,u,®)as  its  coordinate  in 

A 

(ju,u,co)  system.  Then  (7 ,  the  directional  vector  of  the  line 
integral  along  SA  can  be  uniquely  determined  by  i>,  v  and  a. 
So  the  cone-beam  projection  along  line  SA  can  be  described  by 
g(<b,v,a) .  Now  Grangeat's  formula  can  be  expressed  as 


op 


SC 


SO 


SO 


SA 


-g(<b,u,co)dq 


where  p  =  ucosa  +  wsina,  q  =  -usina  +  fi;cosa  . 
Assume  that 


(5) 


SO 

G  (i>,v,a)  = - 2~g(i>,u,a>)  ,  (6) 

|sa| 

so  we  get  the  new  equation  from  equation  (6)  as  follows: 

I _ ,i2 

—R(S,  p)  = L  f+”[  cosa — G  (<t>, u,  co)  +  si n  a —  G  (<I>,  u, a) ]dq 
dp  |S0|  dv  da 

(7) 

Since  the  partial  derivatives  —  G[<b,v,a)  and 

dv 


—G(i>,v,a)  need  to  be  calculated  only  once  in  the  whole 

da 

reconstruction  process,  the  computational  complexity  is  greatly 
reduced. 


C .  P  rojections  R  ebinning 

The  rebinning  of  projections  is  to  find  the  transform  functions 
among  the  parameters  on  the  left  side  and  the  right  side  of  equation 
(8).  That  is  for  each  point  ( 0,<p,p )  in  the  Radon  space,  calculating 
the  corresponding  S,  the  focal  point  position;  a  ,  the  direction  of 
the  integration;  and  p  ,  the  distance  from  the  line  of  the  integration 
to  the  origin  of  the  coordinates. 

To  our  volume  scanning  structure,  we  adopt  the  geometry  with 
the  positions  of  detectors  and  x-ray  source  reversed.  In  the  case  of 
image  reconstruction,  the  circle  and  arc,  which  the  detectors  are 


located,  can  be  regarded  as  the  orbit  of  cone- beam  vertexes;  and  the 
trough-like  x-ray  source  can  be  treated  as  detector  pixels  on  the 
trough-like  surface.  The  reconstruction  algorithm  of  orthogonal 
circle-and-arc  orbit  volume  CT  can  be  utilized  in  our  system.  Hence, 
we  adopt  the  backward  mapping  method  for  rebinning  projections 
here  [7], 

We  assume  that  any  Radon  plane  intersects  the  circle  orbit  at 
two  points,  named  S3  and  S2,  and  the  position  arrangement  for 
S2— 0  is  counterclockwise  if  one  observes  from  above  the  circle 

orbit.  Second,  let's  suppose  that  the  angle  between  OS1  and  the  x 


axis  is  </>cl  and  that  between  OS2  and  the  X  axis  is  ^c2.Then, 
for  a  given  point  [6,(p,p)  in  the  Radon  space,  </>cl  and  ^c2  can  be 

calculated  directly  from  the  coordinates  of  S2  and  S2,  respectively. 

(1)  Rebinning  from  the  circular  orbit: 

For  a  given  [0,tp,p),  p  and  a  can  be  solved  exactly  for 
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(2)  Rebinning  from  the  arc  orbit: 
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(8  b) 
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(10  b) 


Next  step,  we  need  to  find  the  mapping  of  the  projection  data 
from  the  plane  to  the  trough-like  surface.  In  Figure  8,  for  a  point 
(&,a,p)  in  the  detector  plane  there  is  a  point  (®,/?,t)  in  the 
trough-like  surface,  where  /?  and  t  can  be  calculated  via  the 
following  equations: 

r,2  +r,2  -x2 

B  =  n  -  arccos^ - - - 

2r2r3 

xplcosal 

t=  ,  1  '^= 

y(r2  +r3)2  +  p2  sin2  a 

where  x  is  the  positive  solution  of  the  following  equation: 

r22  =  x2  +  rj  -  2 xr3  cos (arctg  1111^)  (12) 

r2r3 

So  far,  based  on  the  first  derivative  of  the  Radon  transform 


(11  a) 
(11  b) 


—  R(S,p )  from  the  rebinning  process,  calculation  of  the  second 

dp 

derivative  can  be  accomplished  by  convoluting  —R(S,p)  with  a 

dp 

1-D  derivative  filter.  The  object  f  (x)  can  be  reconstructed  by  back 
projection,  as  indicated  in  equation  (4). 

D.  3-D  Reconstruction  Algorithm 

Once  the  cone- beam  projections  are  obtained,  the  procedure  to 
reconstruct  the  3D  images  from  the  projections  is  descried  as 
follows: 

1)  Pre-weighting  cone-beam  projections,  using  equation  (6); 

2)  Calculating  the  partial  derivatives  of  the  pre-weighted 
projections,  using  equation  (7); 

3)  Rebining  the  first  derivative  of  Radon  transform,  using 
equations  (8),  (9)  and  (10); 

4)  Calculating  the  second  derivative  of  Radon  transform  , 
using  equation  (11)  and  (12); 

5)  Reconstructing  the  images  by  inverse  Radon  transforms, 
using  equation  (4). 

IV.  Simulation  experiment 

We  make  simulation  experiment  with  Shepp-Logan 
phantom  to  verify  our  3-D  imaging  system  and  its 
reconstruction  algorithm,  in  the  experiments,  we  place  point 
detectors  3°  apart  on  the  circle-plus-arc  orbit.  The  total 
number  of  detectors  is  130.  This  is  equivalent  to  cone-beam 
x-ray  scanning  at  130  positions  and  gaining  130  groups  of 
projection  data.  A  transverse  slice  of  the  Shepp-Logan 
phantom  at  z=-38.5  mm  is  in  Figure  5(a)  for  reference,  which 
is  directly  calculated  from  the  phantom  data.  Figure  5(b)  is 
the  same  slice  reconstructed  from  projection  data.  Their 
corresponding  profiles  are  shown  in  Figure  5  (c)  and  (d).  The 
artifact  in  the  reconstruction  image  are  due  to:  (1)  sampling 
rate  (the  number  of  detectors  on  the  circle-plus-arc  orbit),  (2) 
the  error  from  two  phases  of  rebinning  of  projection  data  to 
Radon  space. 

V.  Conclusion 

The  volume  CT  with  a  trough-like  source  and  point-detectors 
proposed  in  this  paper  has  some  innovations  as  follows: 

(1)  A  novel  geometry  of  raster  scanning  x-ray  source  that 
simplifies  machine  structure. 

(2)  Using  certain  distributing  detectors,  the  imaging  system 
could  gain  full  projections  through  the  trough-like  x-ray  source 
rotation  only  a  single  time.  So  it  achieves  optimal  combination  of 
the  following  aspects:  minimize  rotating  times,  Compton  scatter  to 
be  restrained,  and  image  reconstruction  speed  to  be  improved. 

And  from  above  simulation  experiments,  we  can  see  that  the 
direct  volume  CT  scanning  structure  proposed  in  this  paper  achieves 
not  only  faster  scanning  speed  but  also  more  precise  reconstruction 
image  than  the  methods  in  the  literature.  Moreover,  this  structure 
does  not  need  to  move  the  patient  during  scanning,  which  is  suitable 
for  image-guided  surgery.  Besides  that,  this  system  specially  suits 
for  micro-animal  imaging,  which  can  be  used  to  check  the  micro 
focus  of  infection,  track  the  curative  effect  of  new  medicine  or  gene, 
and  so  on. 


(a)  (b) 


(c)  (d) 

Figure  5.  The  comparison  of  original  image  with  reconstruction 
image  at  certain  position.  Figure  (a)  and  (b)  are  the  slice  images 
of  Shepp-Logan  module  and  reconstruction  image  at  z=80. 
Figure  (c)  and  (d)  are  the  profiles  at  the  position  indicated  by 
dashed  lines  in  figure  (a)  and  (b). 
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